Electron-phonon coupling in DyFeO,3 revealed by infrared spectroscopy 
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We have investigated crystal field and phonon dynamics of the multiferroic orthoferrite DyFe03 
via polarized infrared spectroscopy. Reflectance of single crystals was measured in the far- to mid- 
infrared spectral range at range of temperatures from 10-295 K. We observe a strongly anisotropic 
phonon spectrum which differs from earlier lattice dynamical calculations in its symmetry, as well 
as a mode with significant and unusual temperature dependence that we interpret as a coupled 
phonon-crystal-field excitation. 



Multiferroic materials have long been of interest due to 
the rich variety of magnetically- and electrically-ordered 
phases they exhibit, and serve as excellent model systems 
for a range of spin-charge coupling phenomena. Further 
interest arises from the prospect of device applications 
which exploit magnetoelectric coupling to exert electric 
control over magnetic states, and vice-versa. Realiza- 
tion of such applications depends firstly on the energy 
scales associated with the ordered states; they must be 
large enough for macroscopic polarization effects to be 
robust at operational temperatures. Of the materials 
demonstrated to exhibit this quality, BiFe03 has been 
the most studied. Rs ordering temperatures T/v ee z = 640 
K [T] and Tc U ne — 1090 K [5] are remarkably high with 
respect to room temperature. The other criterion for 
device functionality, however, is strong magnetoelectric 
coupling, lattice-mediated interactions between the spin 
and charge polarizations which permit their mutual con- 
trol [5]. This coupling is found in magnetically induced 
multiferroics such as CuO [I]. Since BiFeC>3, an im- 
proper ferroelectric, derives its polarization from a lone- 
pair dipole phenomenon, and not from distortions in the 
ionic lattice, it does not possess the necessary coupling 
for applications, leading to a search for new materials 
that exhibit stronger magnetoelectric coupling. 

Recently, attention has fallen to the rare-earth or- 
thoferrite (RFeC>3) series [6], whose orthorhomically- 
distorted perovskite structure contains two magnetic sub- 
systems and tilted oxygen octohedra, providiving am- 
ple opportunities for magnetoelectric coupling via the 
Dzayloshinskii-Moriya interaction [7 . SmFeOa was re- 
cently shown to be an improper ferroelectric, with ME 
coupling over a range of temperatures [5]. DyFeOs has 
been demonstrated [5] to exhibit ME coupling, but only 
at temperatures below the rare-earth Neel temperature 
T^ = 4 K, and in the presence of an applied magnetic 
field. A spin reorientation transition at 37 K, however, 
suggests that higher energy scales may be at play. 

Given the importance of the lattice in mediating in- 



teractions between the magnetic and electric polariza- 
tions, a detailed understanding of the lattice vibrational 
modes and their couplings to spin and charge symme- 
tries is essential to the advancement of multifcrroic-based 
technologies. Infrared spectroscopy is an especially effec- 
tive tool for exploring the lattice spectra; it offers high 
energy resolution and sensitivity for detecting infrared- 
active lattice vibrations, and reveals interactions between 
phonons and other dipolar excitations, both electric and 
magnetic. 

Here we report measurements of the reflectivity of 
DyFeOs in the infared frequency range, and utilize an 
multi-oscillator model fit to parameterize the observed 
absorptions. We find that the experimental eigenfrequen- 
cies and symmetries are not in agreement with those cal- 
culated in previous work on orthoferrites [9 . For elec- 
tric fields along the orthorhombic [1 1 0] crystal axis we 
observe the low-temperature emergence of an absorption 
corresponding to a crystal field level transition; this mode 
is coupled to the lowest frequency phonon, leading to 
notable and highly unusual shifts of eigenfrequency and 
linewidth. 

The single crystals of DyFeOs measured in this study 
were grown by a flux method by J. P. Remeika at Bell 
Laboratories. A 2x2x2.5 mm 3 sample was selected for 
magnetization measurements, while a much larger crys- 
tal from the same batch, with a mirror-like, flat natural 
face, was used for optical and x-ray experiments. X-ray 
diffraction confirms the orthorhombic structure, with lat- 
tice constants a =5.3 A, 6=5.6 A, and c=7.6 A. Sharp 
diffraction peaks indicate that the sample is single-phase, 
without garnet inclusions or flux contamination. Magne- 
tization was measured in a SQUID magnetometer (Quan- 
tum Design MPMS) with magnetic fields along the the 
principle pseudo-cubic axes. 

Infrared reflectance spectra were recorded at near- 
normal incidence for frequencies oj — 50-5000 enr 1 and 
temperatures T — 10-300 K. Incident electric fields were 
polarized along either the [1 1 0] or [1 1 0] axes, corre- 
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sponding to the unit axes of the high-temperature cubic 
unit cell, rather than those of the orthorhombic structure. 
(See inset to Fig. 1.) Since the tilting of oxygen octohe- 
dra is relatively small the crystal retains a pseudocubic 
symmetry. As we show below, these measurement axes 
provide the best basis for understanding the electrodyn- 
mic response. 

Magnetic susceptibility xm is shown in Fig. 1(b) as a 
function of temperature for magnetic field H = 1000 Oc 
applied parallel to the [1 10], [1 10], and [0 1] axes. 
Below T^e the Fe ions are ordered antiferromagnetically 
in the [1 0] (G-type) and [0 1 0] (A-type) directions, 
with a small sublattice canting producing a weak ferro- 
magnetic moment (WFM) in the [0 1] direction ®[H]. 
At the spin reorientation transition temperature Tg R e = 
37 K the spin structure rotates from G x A y F z to A x G y C z , 
and the WFM disappears. Dy ions order at T^ v = 4.5 
K. 

Reflectivity of DyFeC>3 in the far-infrared frequency 
range, shown in Fig. 2 for selected temperatures, is dom- 
inated by lattice vibrational modes. These modes gener- 
ally can be assigned to three different types of motion: Dy 
external modes at the lowest frequencies, Fe-0 bending 
modes in the intermediate range, and oxygen stretching 
modes at higher frequencies. At room temperature the 
differences between spectra taken with [1 1 0] (Fig. la) 
and [1 -1 0] (Fig. lb) polarizations are readily apparent, 
with 16 phonons observed in the former and 12 in the 
latter. As the temperature is lowered to T = 10 K, all 
modes in the [1 1 0] spectra harden and narrow slightly 
due to nonlinear effects. The [1 1 0] spectra, however, are 
strongly renormalized at low temperature. In particular, 
a new mode appears at 140 cm" 1 for T <180 K. 

Spectra recorded at these two polarizations form a lin- 
early independent basis; that is, spectra measured at in- 
termediate polarizer angles contain an admixture of the 
basis spectra. This strong association of the phonon 
spectra to the pseudocubic unit cell is at odds with the 
dynamical lattice calculations of Ref. 9, which predict 
phonon modes of type Bi u , i?2u, and B^ u polarized along 
the orthorhombic crystal axes. This discrepancy makes 
phonon assignment difficult, and highlights a need for 
further effort in calculating the lattice dynamics. 

In order to quantify the electrodynamic response of 
this system we fit the reflectance spectra to a sum of 
Drude-Lorentz oscillators of the form 

N 2 

where uj Ptn is the plasma frequncy, wo,n is the transverse 
resonance frequency, and 7„ is the linewidth of the nth 
oscillator. A fit to the four-parameter product form of the 
Drude-Lorentz oscillator [TO], which takes into account 
LO-TO phonon splitting, was also carried out; however, 
since the fit was not substantially improved, the more 
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FIG. 1: (a) Crystal structure for orthorhombically distorted 
perovskite DyFe03. Orthorhombic axes are rotated 45° rela- 
tive to the cubic unit cell, (b) Magnetic susceptibility M/H, 
collected with H = 1000 G oriented along the principle pseu- 
docubic axes. 

intuitive 3-parameter fit results will be discussed here. 
Best-fit eigenfrequencies and linewidths are shown in Fig. 
3, scaled to their values at 180 K. As seen in Fig. 3a, 
temperaure-induced changes to [1 1 0] phonon modes are 
less than 2%, and are qualitatively similar for all modes 
in the far infrared region. 

For E || [1 1 0], however, the lowest-frequency modes 
change dramatically with temperature. The lowest- 
frequency mode, labeled mode TOl, is only resolved be- 
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FIG. 2: (a) Infrared reflectivity for DyFeOa with E || [1 1 0] 
polarization, (b), (c) Reflectivity for [1 1 0] polarization. 

low T m 180 K, and exhibits an 7% downshift in eigenfre- 
quency upon cooling to 10 K. The next highest mode, in 
contrast, hardens significantly as temperature decreases, 
with a 7% increase in transverse frequency. 

The temperature dependence of phonon damping is 
also anomalous. Unlike the case for [1 1 0], where all 
modes exhibit nearly identical narrowing of 20-50%, for 
[T 1 0] a range of damping values are observed, with 
the greatest changes occurring in the lowest-frequency 
modes. For T02, the low-frequency mode that hard- 
ens as temperature is lowered, 7(10 K) drops as low as 
0.157(180 K). The linewidth of TOl is especially pecu- 
liar, increasing slightly after TOl is first observed, then 
decreasing sharply below T = 135 K. 

In interpreting the origin of the TOl mode observed 
here, recall that there are many electronic processes 
which produce far-infrared absroption signatures. We 
now discuss several of these possibilities and conclude 
that the explanation most consistent with the data is 
that of a crystal field transition coupled to two nearby 
phonons. The first possibility is that TOl is a lattice 
mode. Indeed, it has comparable frequency and oscillator 
strength to nearby phonon modes, but its sudden onset 
and narrowing at the intermediate temperature scale of 
135-180 K would be difficult to understand in the absence 



FIG. 3: Multi-oscillator fit results for E || [T 1 0] polarization. 
Fit parameters (a) eigenfrequency and (b) linewidth are scaled 
to their values at T= 180 K. (c) Dielectric strength for the 
two lowest frequency modes, as well as total dielectric strength 
from all infrared modes. 



of a structural phase transition. Such a transition has not 
been previously reported in related rare-earth orthofer- 
rites; thus, this assignment is highly unlikely. Dramatic 
decreases in phonon transverse frequency and oscillator 
strength are commonly seen in systems with spontaneous 
dielectric polarization, as lattice vibrations with eigen- 
vectors corresponding to the polarization vector become 
unstable. However, a study of the dielectric constant in 
DyFeOa only found spontaneous polarization for T < 4 K 
and magnetic fields greater than 24 kOe [8] . The present 
optical data agree with this result, as the Lydanne- Sachs- 
Teller sum rule dictates that such changes in the dielec- 
tric strength of finite-frequency oscillators must be offset 
by an increase in the low- frequency dielectric value eo- 
As seen in Fig. 3(c), the sum of the dielectric strength 
of all observed oscillators remains constant, and eo is un- 
changed. In fact, most of the dielectric strength gained 
by the new mode TOl is offset by a decrease in the 
strength of T02. 

Given the high ferromagnetic ordering temperature 
and the antiferromagnetic and spin-reorientation transi- 
tions at low temperatures, it is worth considering whether 
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FIG. 4: Crystal field level populations calculated using Boltz- 
mann statistics 



the TOl mode could be magnetic in origin. How- 
ever, the large dielectric strength and relatively high fre- 
quency are inconsistent with an assignment of a magnetic 
dipole transition. Electrically-excited magnons (electro- 
magnons) have been observed in the related perovksitc 
manganite (RMnOa) compounds; these excitations are 
most commonly observed at sub-THz frequencies and in 
regions of the temperature-magnetic field phase diagram 
with magnetically-induced ferroelectricity, but they have 
been found in non-multiferroic phases, and at higher fre- 
quencies and temperatures OLH]- In the latter cases the 
excitations are often related to two-magnon processes, 
emerging out of spin fluctuations at temperatures sever- 
als times greater than the ordering temperature [HI H2] . 
In the case of DyFeC>3, however, the spin structure lacks 
noncollincarity, and thus does not produce the dielectric 
polarization and spin-lattice coupling central to an elec- 
tromagnon. Also, the intermediate temperature scale at 
which TOl turns on (180 K) is commensurate with nei- 
ther the rare-earth T Ned (4.5 K) nor the Fe T c (745 K), 
and the mode is insensitive to the change in magnetic 
symmetry at T$r — 37, further discounting the possibil- 
ity of an origin in a collective magnetic phenomenon. 

We propose that the TOl mode is the signature of an 
electric dipole transition between crystal field-split levels 
of the electronic groundstate. Earlier optical studies of 
DyFeOs [14] have revealed 6 -ffi5/2 — s- 6 P5/2 transitions 
with crystal field levels I, II, and II at 52, 147, and 225 
cm -1 , respectively. The second level (147 cm -1 ) is close 
in energy to the TOl mode, and matches its extrapo- 
lated value at 200 K (near the temperature where TOl 
is first resolved). These levels are thermally populated 
according to Boltzmann statistics, and their relative pop- 
ulations Pi are have been calculated and plotted in Fig. 
4 as a function of temperature. Transitions from the 
ground state to level i begin when Pq — Pi > 0; this onset 



is calculated to be 125 K for level II, in rough agreement 
with the data. 

We suggest that the anomalous temperature depen- 
dence and of both the TOl and T02 modes, as well as 
the transfer of spectral weight between them, represents 
the realization of a coupled crystal-field-phonon excita- 
tion 18-20J. Such a coupling can only be observed when 
a rare-earth crystal field (CF) transition lies close in en- 
ergy to a phonon of the same symmetry, and the two 
resultant excitations of the composite system assume a 
mode-repulsion character, diverging in energy with in- 
creasing interaction strength. CF-phonon coupling is 
made allowed in this system due to the distorted crys- 
talline electric field at the orthorhombic Fe sites. De- 
tailed calculations [T8H2Tj] based on a simple two-level 
model have shown that, for coupled two-peak structures 
in Raman scattering data, the unrenormalized mode po- 
sitions and intensity ratio can be used to determine the 
bare CF excitation and phonon frequencies ujcf and u p h, 
as well as their interaction strength V. Here we extend 
this analysis by analogy to modes of infrared activity to 
extract locf, ^>ph, and V via 
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V = |w ,i - U>0,2\: 
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where R — ojp^/ui^. The bare values have been cal- 
culated from parameters obtained from the Lorentzian 
oscillator fits to the reflectivity and plotted in Fig. 5. 
The interaction strength V increases as the tempera- 
ture is lowered, reflecting the greater splitting between 
modes and relative gain of spectral weight by the TOl 
mode. The bare frequencies, in contrast to the ob- 
served ones, are nearly constant with temperature, with 
ujcf = 147— 149 enr 1 closely matching the value of 147 
cm -1 from earlier optical studies [T4"] . 

In summary, we have observed phonon-crystal-field 
coupling in the muliferroic perovskite DyFeOs via in- 
frared spectroscopy. Calculation of the bare CF excita- 
tion and phonon energies yields excellent agreement with 
earlier optical meaurements, adding credence to this in- 
terpretation of the infrared data. These findings yield 
new insight into the nature of the interaction of the crys- 
tal lattice with electronic degrees of freedom, an essential 
ingredient in the realization of new multiferroic phenom- 
ena. More detailed lattice dynamical calculations will im- 
prove our understanding of the symmetries of the modes 
involved and their relation to the electronic and magnetic 
structure. 
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FIG. 5: (a) Bare crystal field and phonon frequencies 
extracted from experimentally observed frequencies. (b) 
Crystal-field-phonon coupling parameter V. 
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